Dymamical Casimir emission from polariton condensates 



Selma Koghee 1 and Michiel Wouters 1 

1 TQC, Universiteit Antwerpen, Universiteitsplein 1, B-2610 Antwerpen, Belgium 

(Dated: March 6, 2013) 

We study theoretically the dynamical Casimir effect in an exciton-polariton condensate that 
is suddenly created by an ultrashort laser pulse at normal incidence. As a consequence of the 
abrupt change of the quantum vacuum, Bogoliubov excitations are generated. The subsequent 
evolution, governed by polariton interactions and losses, is studied within a linearized truncted 
Wigner approximation. We focus in particular on the momentum distribution and spatial coherence. 
The limiting behavior at large and small momenta is determined analytically. A simple scaling 
relation for the final condensate depletion as a function of the system parameters is found and the 
correlation length is shown to depend linearly on the condensate depletion. 



Interacting systems possess, as a profound consequence 
of quantum mechanics, a rich ground state structure. An 
immediate implication is the dynamical Casimir effect: 
when some parameter is abruptly varied, the system no 
longer finds itself in the new ground state, causing the 
appearance of excitations. Since the original idea of gen- 
erating photons by accelerating two mirrors with respect 
to each other [TJ [2], the dynamical Casimir effect has 
been studied in contexts as diverse as particle production 
in a time dependent cosmological background [3-5], su- 
perconducting circuits [6j [7] and ultracold atomic gases 
[8HTU]. In the latter system, the nontrivial vacuum is 
that of the (weakly) interacting Bose gas, which features 
quantum depletion. 

A different realization of the weakly interacting Bose 
gas is a microcavity containing exciton polaritons [TT] . 
Their advantage over atomic gases is that exciton polari- 
tons have a photonic component, offering great flexibil- 
ity for optically creating them in a specific state. This 
feature was e.g. exploited in the experiments on polari- 
ton superfluidity, where the polariton velocity could be 
varied by simply changing the angle of incidence of the 
excitation laser [T2] . 

This advantage has recently led to several proposals 
to utilise exciton polaritons for the experimental obser- 
vation of Hawking radiation p~3j [14] , a phenomenon in- 
timately related to the dynamical Casimir effect. In this 
Letter, we analyze a simpler manifestation, directly re- 
lated to the quantum depletion of the interacting Bose 
gas. We propose to create polaritons by a pulsed laser at 
normal incidence, which results in a state with all parti- 
cles at zero momentum. Since this is not the true ground 
state of the many body system, quasi-particles are ex- 
pected to appear. In other terms, by the sudden cre- 
ation of a coherent polariton population, the vacuum is 
changed from the trivial one to the Bogoliubov vacuum, 
leading to dynamical Casimir emission. The change of 
the dispersion due to the presence of a condensate was 
probed in the experiments of Kohnle et al. [15 by a four 
wave mixing pump probe experiment. We suggest here 
to study the dynamics without the probe pulse. 

An important aspect of polariton physics in actual ex- 



periments is the finite polariton life time. This brings 
the system in an inherently out of equilibrium situation 
and limits the possibilities to observe quantum effects. 
On the other hand, it offers a unique laboratory to study 
the interplay and competition between quantum correla- 
tions and dissipation, that may teach us lessons that are 
relevant for cosmological studies [T6] . A basic question 
is how the achievable quantum depletion depends on the 
system parameters. We will show that within the Bogoli- 
ubov approximation, the dynamical quantum depletion 
of a one dimensional polariton wire [17] scales linearly 
with a "blockade" parameter U/j [18] , where 7 is the 
line width and U is the interaction energy of two po- 
laritons that are localized within the condensate healing 
length. 

The coupling between the polariton fluid and the out- 
side world leads both to particle losses and to the en- 
trance of quantum fluctuations into the system. This 
coupling can be described accurately with the master 
equation. For a quadratic Hamiltonian, an equivalent 
formulation in terms of the Wigner distribution of the 
polariton field can be constructed. This approach comes 
down to mapping the quantum master equation to a set 
of classical stochastic differential equations to sample the 
Wigner distribution function. The stochastic differen- 
tial equation for the polariton field reads in the so-called 
truncated Wigner approximation [19] 
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where e(— iV) gives the kinetic energy of a free particle, 
g is the strength of the contact interaction, and AV is 
the volume of the grid cell necessary in the derivation 
of the truncated Wigner approximation. The stochas- 
tic term dW is white Gaussian noise with amplitude 
variance (dW*(x, t)dW(x.', t)) = 2dtS yi ^ and a random 
phase. Expectation values with respect to the Wigner 
distribution of the stochastic fields correspond to the 
symmetric operator averages, e.g. (^(x, ^^(x', t') + 
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Bogoliubov matrix 



Figure 1. When a short laser pulse (1) resonantly creates a 
polar it on condensate (2), a sudden change of trivial vacuum 
to the Bogoliubov one, results in the generation of dynamical 
Casimir emission (3) of Bogoliubov excitations. 



^(x'^'^CM)) = (0*(x,£)0(x',£)) w . This allows for 
a direct computation of all equal time correlation func- 
tions [20] . 

The action of a resonant excitation pulse at k = can 
be included in Eq.Q by adding a pumping term Fl(x, t) 
on the right hand side. However, if the pulse is very short, 
it can be replaced by a shift of the initial condition of the 
stochastic variable 
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where 77 (x) is the normally distributed uncorrelated vac- 
uum noise with variance (^*(x)^(x / ))v^ = 1/2 £ x?x /. The 
subsequent dynamics is no longer influenced by the ex- 
citing laser, which merely fixes the initial condition. This 
situation is very different from both continuous wave res- 
onant and nonresonant excitation schemes, where an in- 
terplay between driving and dissipation takes place at 
all times. Compared to the continuous wave excitation 
regime, the dynamics of the polar iton gas in the present 
scheme resembles the standard Bose gas more closely. 

Only in the case of a quadratic Hamiltonian, the map- 
ping between the quantum and classical problem becomes 
exact. Since we will restrict our calculations to the Bo- 
goliubov approximation, the Wigner forumulation does 
not introduce any further errors. Let us therefore de- 
compose the polariton field in a condensate (k = 0) and 
fluctuation component; </>(x,£) = (j) c {t) + £(/>(x, t). In the 
equation of motion for the condensate, the effect of the 
interactions with the noncondensed particles and the ef- 
fect of fluctuations are negligible within the Bogoliubov 
approximation. Thus, the time evolution of the conden- 
sate component is given by [15] 



c(t) = V n c(t) exp 
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where the condensate density n c (t) = n^e~ 7t 
exponentially and its phase is given by the time inte- 
grated blue shift c (t) = g J Q n(t f )dt r . Up to first order, 
the motion of the fluctuations is governed by the familiar 
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where $ k = [e^^d^k) e^£0*(-k)] T and dS k (t) = 
[dW(k,t) - dW*(-k,£)] T , with the correlation of 
the noise in momentum space (dW*(k, t) dW(k', t)} = 
2dt <5 k;k /. 

Due to the time dependence of the Bogoliubov ma- 
trix, an analytic solution for the fluctuations <£ k is not 
straightforward to compute, but for small and large wave 
vectors, limiting expressions can be found analytically. 
Let us start with the limit k — 0. In this case, the 
time dependence of the real part of the Bogoliubov ma- 
trix comes down to a rescaling with a factor exp(— 7t/ft), 
whereas the imaginary part is proportional to the iden- 
tity matrix. Hence, the Bogoliubov matrices at different 
times commute with each other. The differential equa- 
tions are then easily solved, yielding for the momentum 
distribution 
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where we have used the relation between the momentum 
distribution and the expectation value of the stochastic 
fields that sample the Wigner distribution W. At large 
times, Eq. (p| reduces to the simple scaling n(k — >• 0, t — >• 
00) = 2(gn^/j) 2 exp(-jt/h). 

For the other limit of large fc, we can use the sud- 
den approximation from [21] to compute the number of 
Bogoliubov excitations that are created. At zero temper- 
ature, this approximation yields UB^k) = \v(k)\ 2 , where 
v(k) = \/[e(k) + gnW]/[2hu B (k)] - 1/2 is the usual Bo- 
goliubov v coefficient which depends on the dispersion as 
huj B (k) = y/e(k) [e(k) + 2griW\. The polariton momen- 
tum distribution is related to the quasi-particle occupa- 
tion as 

(^(k)V'(k)) = u(k) 2 (b\k) 6(k)) + v(k) 2 (b(-k) 6 f (-k)) 

+u(k)v(k)(b(k) 6(-k)) + u(k)v(k)(tf (-k) 6 f (k)), 

(7) 

The terms on the second line represent an oscillating 
contribution, that we will not take into account in the 
following. Using the initial occupation of the Bogoli- 
ubov excitations, and discarding losses one obtains for 
the momentum distribution n(k, t) = 2\u(k)\ 2 \v(k)\ 2 . In 
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the crudest approximation, losses can be included by an 
exponential decay: 
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The wave number k* where the cross-over between the 
small (|6| and large Q momentum behavior takes place 
can be found by equating both expressions. In the more 
interesting limit where gn^ ^> 7, this yields 
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In the limit where the decay rate becomes much smaller 
than the interaction strength, this wave vector tends to 
zero and the expression ^ becomes accurate down to 
very small wave vectors. 

For the behavior of the momentum distribution around 
k*(t), we have not found an analytic formula, but the Bo- 
goliubov dynamics ([5| was solved numerically instead. 
To this purpose, we have introduced the Green's func- 
tion of the time dependent linear system G^(t : t f ). The 
momentum distribution can then be expressed as 
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The Green's function itself was computed by discretising 
time interval [t = t/v,^ = t\] in small steps At 
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The momentum distributions obtained with this nu- 
merical scheme are presented in Fig. [2] We have rescaled 
them to the condensate density in order to cancel the 
trivial exponential decay. The red dotted line repre- 
sents the expression (|8|, that is seen to match very well 
for sufficiently large momenta, justifying the sudden ap- 
proximation in this regime. For the smaller interaction 
strength (a), the momentum distribution decays as /c -4 , 
where for the larger interaction strength (b), its decay 
is dominated by the intermediate k~ 2 region. The black 
solid line represents the prediction §6§ for t —> 00, which 
matches perfectly with the numerical result at late times. 

From the full momentum distribution, we can now 
compute the total number of particles that are produced 
by the dynamical Casimir effect. This result will de- 
pend on the dimensionality of the system. Here, we will 
evaluate the particle production for a one dimensional 
polar iton wire, where quantum fluctuations are most im- 
portant due to the high density of states at low momenta. 
At large times, the ratio of particles in nonzero momen- 
tum states to the condensate density n c (t) can be written 
as 
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Figure 2. The momentum distribution n(k) of the polari- 
ton gas due to the dynamical Casimir effect. The full lines 
show the numerical result for various times (t=0. 2,0. 5, 1,3,8). 
The red dotted line represents the asymptotic result (J8|; the 
thin black line is the prediction ^ for small k at large 
times. We have chosen g — 0.01 /xm meV, 7 — 0.05 meV, L — 
100 /xm, Ti — 1, m - 
panel: = 10. 



1. Upper panel: gn^ /j = 1. Lower 




Figure 3. The total quantum depletion as a function of the 
initial blueshift. The circles show the numerical result; the full 
line is the analytical prediction. We have taken an interaction 
strength g = 0.01/xmmeV and a decay rate 7 = 0.05 meV. 



where we have used (|6| for the number of particles cre- 
ated in the small momentum states and we have taken 
the continuum limit. By approximating the integrand in 



Eq. (12) by 1 for k < k* and by using the expression (|8j) 
for k > fc*, we find that the particle production scales in 
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the limit gn^ > 7 as 
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where the constant can be computed numerically from 
the full momentum distribution: C « 1.5. At large times, 
the depletion (13) can be written as Sn/n « 0.77 g/(^j), 



where £ = h/ ymgrJ^ is the healing length of the con- 
densate immediately after its creation. This shows that 
the particle production is proportional to the "blockade" 
parameter C//7 [H] , where U = g/£ is the interaction 
energy of two polaritons localized within the condensate 
healing length. Typical attainable values of the heal- 
ing length are one micron, so that with an interaction 
strength g = 0.01meV/im, a linewidth of 0.05 meV is 
sufficient to obtain more than ten percent of quantum de- 
pletion. Figure |3] compares the analytical estimate (13) 
with the depletion computed from the full numerical mo- 
mentum distribution. It turns out to be already accurate 
for moderate values of the parameter gn^/j > 1. Al- 
though large values of this parameter seem attainable 
with state of the art microcavities [22 , this regime goes 
beyond the scope of the present work. When gn^ /j 
becomes too large, interactions between excitations can 
no longer be neglected, so the Bogoliubov approximation 
breaks down. 

Experimentally, the quantum depletion is directly ac- 
cessible through a measurement of the first order spa- 
tial coherence g^ 1 \x — x') = (ip^(x, t)ip(x', t))/n c (t). The 
fluctuations that are created because of parametric scat- 
tering lead to a drop in the first order coherence function 
equal to the quantum depletion: g^\x — » 00) —g^\x —> 
0) = Sn/n c . We restrict again to a one dimensional wire, 
as we did for Sn/n c . 

The numerical results for the first order coherence are 
presented in Fig. [4j Since the gas is assumed to be 
created by a perfectly coherent laser, no fluctuations 
are present initially. The correlation length can be es- 
timated from the relation £ c (t) oc l/k*(t). The red cir- 
cles in Fig. [4] show this prediction to work well. From 
this relation, it is seen that the coherence length scales 
at short times as £ c (t) ex 2t\J gn(°) Jm. This is indeed 
what one expects for fluctuations that spread out at the 
speed of sound c = \J gn c (t)/m. At long times, we obtain 
t c (t — >• 00) = 7.8 \/ gnW Jm ft/7, where the constant 7.8 
was determined numerically. This means that the effec- 
tive propagation time of the fluctuations is the life time 
r = 7 _1 ; their average speed is proportional to the initial 
sound velocity. Finally, it is interesting to note in Fig. [4] 
the linear relation between the time dependent depletion 
and the coherence length: £ c (t) oc Sn(t)/n c (t). By com- 
bining Eqns. ([9| and (13), this relation is readily verified 
analytically. 

The dynamical Casimir generation of Bogoliubov exci- 
tations that we have studied here, is only the first process 
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Figure 4. The spatial coherence g^\x) as a function of time 
(t — 0.1,0.5, 1,3,8), for the parameters as the corresponding 
panels in Fig. [2] The red circles represent the analytical 
estimate of the correlation length £ c (t) — 3.9//c*(£). The red 
dotted line shows the linear relation between the correlation 
length and the depletion. 



on the route of the Bose gas to thermal equilibrium. For 
long polariton life times, a second stage with collisions be- 
tween Bogoliubov excitations should set in. This opens 
the prospect to explore the rich physics of the thermal- 
ization of quantum systems [23] in a new setting. 
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